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ABSTRACT 


This report covers research conducted during the first year of the two- 
year grant. The research, entitled “Application of Lanczos Vectors to Control 
Design of Flexible Structures” concerns various ways to obtain reduced-order 
; mathematical models for use in dynamic response analyses and in control design 
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studies. This report summarizes research described in the following reports and 
papers that were written under this contract. 

• Su, Tzu-Jeng, A Decentralized Linear Quadratic Control Design Method 
for Flexible Structures, Ref. [5] 

• Su, Tzu-Jeng, and Craig, Roy R. Jr., “Controller Reduction by Preserving 
Impulse Response Energy,” Ref. [7] 

• Su, Tzu-Jeng, and Craig, Roy R. Jr., “Substructural Controller Synthe- 
sis,” Ref. [8] 

• Su, Tzu-Jeng, and Craig, Roy R. Jr., “Substructuring Decomposition and 
Controller Synthesis,” Ref. [9] 

• Craig, Roy R. Jr., and Su, Tzu-Jeng, “A Review of Model Reduction 
Methods for Structural Control Design,” Ref. [10] 

• Su, Tzu-Jeng, and Craig, Roy R. Jr., A Decentralized Linear Quadratic 
Control Design Method for Flexible Structures, Ref. [11] 

• Craig, Roy R. Jr., “Recent Literature on Structural Modeling, Identifi- 
cation, and Analysis,” Ref. [12] 
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Chapter 1 


INTRODUCTION 


This report summarizes research accomplished during the first year of a 
two-year grant on the topic of the application of Krylov and Lanczos vectors 
to the control of flexible structures. Under a previous contract with NASA- 
JSC, entitled “Application of Attachment Modes in the Control of Large Space 
Structures,” [1] Lanczos vectors, and the closely-related Krylov vectors, were 
first applied to problems related to the control of flexible structures. Results of 
this research were published in Refs. [2-4], where the use of Krylov vectors as 
basis vectors and the concept of parameter-matching were combined together to 
develop Krylov model-reduction algorithms. Algorithms for generating Krylov 
vectors for general linear systems and for undamped and damped structural 
systems are presented in these references. 

The design of a controller based on a reduced-order system model can 
lead to three types of control energy spillover: control spillover, observation 
spillover, and dynamic spillover. The combined effect of the three types of 
spillover usually degrades the performance of the controller, when it is applied 
to the full order system, and may even destabilize the closed-loop system. In 
Refs. [2-4] it is shown that, if model reduction and controller design are based 
on a Krylov model, then the control and observation spillover terms can be 
eliminated while leaving only the dynamic spillover to be considered. 
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The goals stated for research under the present grant are: 

1. To develop a theory of reduced-order modeling of general linear systems 
based on the use of Lanczos vectors, and to apply the theory to the 
modeling of flexible structures. 

2. To address numerical issues that arise in the application of Lanczos vec- 
tors to reduced-order modeling, e.g., sensitivity to choice of starting vec- 
tors, loss of orthogonality, etc. 

3. To develop control system design techniques employing Lanczos modeling 
of the controlled and residual systems, considering relevant issues such 
as stability of the closed-loop system, spillover, robustness, and compu- 
tational requirements. 

4. To apply Lanczos-based control system design to general typical prob- 
lems, e.g., optimal co-located velocity feedback, dynamic output feed- 
back, optimal control of finite-time slewing of a beam, etc. 

Much of the research on topics listed under the first two goals above was 
actually conducted just prior to the start of the present grant and is summarized 
in Ref. [1]. A theory of reduced-order modeling of general linear systems and 
of damped and undamped structures (Goal 1) was developed in Ref. [2-5]. 
Steps that may be taken to address several numerical issues that arise in the 
application of Lanczos vectors to reduced-order modeling (Goal 2) are described 
in Ref. [6]. 
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Research that addresses Goals 3 and 4, the development of control system 
design techniques employing Lanczos vectors, and the application of Lanczos- 
based control system design, is described in Refs. [3,5,7-11]. Three major 
topics were studied: the development of a controller reduction method based 
on Krylov vectors [7], substructuring decomposition and controller synthesis 
[8,9], and a comparison of Krylov-based model reduction with other model- 
reduction techniques [10]. This research is fully described in Refs. [5,11], and 
abstracts of the major published papers are included in Chapter 2 of this report. 

The preparation of a survey chapter entitled “Recent Literature on Struc- 
tural Modeling, Identification and Analysis” [12] was supported, in part, by the 
present grant. 

Finally, Chapter 3 outlines several topics to be studied during the second 
year of this project. 



Chapter 2 


ABSTRACTS OF TECHNICAL PAPERS 


2.1 Controller Reduction by Preserving Impulse Re- 
sponse Energy (Ref. 7) 


This paper presents a controller reduction method that is based on a 
projection subspace called a Krylov subspace. The Krylov subspace is generated 
by a Krylov recurrence procedure. The reduced-order controller is called an 
Equivalent Impluse Response Energy Controller (EIREC) because it has the 
same impulse response energy as the full-order controller. 


The plant to be controlled is described by the standard first-order state- 


space form 

z — Az T Bu + Nzu 

y = Cz + v 

The full-order controller to be reduced is described by 

q = Eq + Fy 

u = Gq 


( 2 . 1 ) 


( 2 . 2 ) 


The impulse response energy is defined as the L 2 energy norm of the impulse 


response of the controller 


f OO 

£ = \\H\\l> = tr[ J H T Hdt ] (2.3) 

where H = Ge Et F is the impulse response of the controller. It is shown that 
the impulse response energy can further be expressed as 

£ = tr[J™ F T e ETt G T Ge Et Fdt] = tr[F T W 0 F ] (2.4) 
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Or, equivalently, 

£ = <r[jT Ge Et FF T e ETt G T dt ] = tr[GW c G T } (2.5) 

where W 0 and W c are the observability and controllability grammians. 


Two algorithms for generating projection subspaces for controller sys- 
tem transformation are presented. One algorithm generates a subspace which 
normalizes the controllability grammian, and the other normalizes the observ- 
ability grammian. The controller is transformed to new coordinates called 
normalized grammian coordinates by using the projection subspaces. Con- 
troller reduction is based upon the representation in the new coordinates. The 
subspace-generating algorithm is a recursive process with either the F matrix 
or the G T matrix as the starting block of vectors. The recursive process is a 


Krylov recurrence procedure which can be summarized by 

Choose L\ 

Li + 1 = M Li 


(2.6) 


The first algorithm uses M = E~ l and Li = F and normalizes the 
vectors L , with respect to the observability grammian W Q . The projection 
subspace is formed as L = [L\ L 2 ...], which satisfies L T W 0 L = /. The L 
subspace is then partitioned into two subspaces L = [L R L T ] with subspripts R 
and T denoting the retained portion and the truncated portions, respectively. 

The reduced-order controller is described by 

'q R = E R q R + F R y 
u = G R q R 

where q = L R q, and where the reduced controller system matrices are given by 
E r = L t r W 0 EL r , F r = L t r W 0 F, and G R = GL R . The F R matrix has nonzero 
elements only in the first block. 
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The second algorithm uses M = E~ T and Li = G T and normalizes the 
vectors L, with respect to the controllability grammian. The reduced-order 
controller system matrices are given by E R = L^EW C L R , F R = L^F, and 
G r = GW c L r . The G r matrix has nonzero elements only in the first block. 

Both algorithms are useful for constructing reduced-order controllers. 
The preference depends on the number of actuators and the number of sensors. 
The algorithms proposed and the reduced-order controllers obtained have the 
following useful properties. 

Property 1: The subspace L is both controllable and observable. If the algo- 
rithm terminates before n vectors are generated, then the full-order controller 
is not minimal. A minimal optimal controller can be produced by projecting 
the full-order controller onto the L subspace. 

Property 2: The reduced-order controller is asymptotically stable if the full- 
order controller is controllable or observable. 

Property 3: If the reduced-order controller is asymptotically stable, then it 
has the same impulse response energy as the full-order controller. 

Property 4: The reduced-order controller matches a set of system param- 
eters. This set of parameters includes the so called low-frequency moments 
and low-frequency power moments. The low-frequency moments are defined 
as GE~'F, i = 1, 2, — The low-frequency power moments are defined as 
F T (E T )- i W 0 E~ j F, and GE~ i W c {E T )-jG T , i,j = 1, 2, .... 

Two examples drawn from other controller reduction literature are used 
to test the proposed controller reduction method. The Equivalent Impulse Re- 
sponse Reduction Algorithm can produce closed-loop system designs with good 



performance. Computationally, the proposed method is economical compared 
with the other controller reduction methods. 


2.2 Substructuring Decomposition and Controller Syn- 
thesis (Refs. 8,9) 

This paper presents a decentralized control design process called Sub- 
structural Controller Synthesis (SCS). A natural decomposition called sub- 
structuring decomposition is used to decompose a flexible structure into several 
substructures. Then, the linear quadratic optimal control design is carried out 
to design a subcontroller for each substructure. The final controller for the 
assembled structure is a global controller, which is assembled from the subcon- 
trollers by using the same assembling scheme as that employed for structure 
matrices. 


Although the method can be applied to structures with more than two 
substructures, for simplicity a two-component structure is used to demonstrate 
the formulation. The equations of motion of the two substructures are repre- 
sented by 




(2.8) 


Mi'ii -f Di&i + K,x, = P iUi 
Vi = ViXi + Wiii 

The dynamics of the the assembled structure (the structure as a whole) is 
described by 


M 'x + Dx -f Kx = Pu 
y = Vx + Wx 

It is shown that there exists a coupling matrix T which relates x a , 
follows: 


( 2 . 9 ) 


xp, to x as 
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and the system matrices of the substructures and assembled structure satisfy 


M = T T 
P = T t 


M a 0 
0 Mp 

Pa 0 

0 Pp 


T,D = T t 


Da 0 
0 Dp 


T,K = T t 


Ka 0 

0 Kp 


V = 


V a 0 
0 V* 


r, 


w = 


w a o 

0 Wp 


( 2 . 11 ) 


For control design purposes, the equations of motion are rewritten in the 
following first-order forms: 


D M 
M 0 


y = [V W] 


x 

x 

X 

X 


-K 

0 


0 

M 


]{:} 


+ 


P 

0 


u 


( 2 . 12 ) 


or 

Sz = Az -f Bu 

y = Cz 

for the assembled structure. Similarly, the equations of motion of substructures 
can be represented by 


(2.13) 


S,z x = AiZi + BiUi 

Vi = C x Zi 


i = a, 0 


(2.14) 


Combination of the two substructure equations above gives the first-order equa- 
tion of motion of the unassembled structure 


Sz = Az + Bu 
y = Cz 


(2.15) 


with the system matrices in the following block diagonal form 


o 

A = 

' Aa O' 

II 

O 

i 

c = 

0 

CJ 

1 

o 

to 

Xu 

? j 1 

o 

XU 

1 

’ L 0 Bp J 


1 

o 

l 



and 
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It can be shown that the system matrices of the unassembled structure and the 
assembled structure satisfy 


S = T t ST A = T t AT B = T t B C = CT 


( 2 . 16 ) 


with 


T = 


T a 0 
0 T a 
Tp 0 
0 Tp j 


The control design method is the LQG theory with the performance 


indices 


1 


Ji = lim -E[xj M,i, + xj M{Xi + uf i = a, fi 


( 2 . 17 ) 


( 2 . 18 ) 


t— ► OO 2 

for substructures being minimized. The subcontroller is in the form 

S t q t = (A + B t G° - F°C t )qi + F t 0 y, 

i = a, p 

Ui = Gfqi 

with F° and G° the optimal control gain and optimal filter gain matrices. In 
a more compact form, the controller equation for the unassembled system is 
written as 

Sd = (A 4 - BG° - F°C)a + F°v 

( 2 . 19 ) 


Sq = (A + BG° - F°C)q + F°y 
u = G°q 


with 


G° = 


1 

CD 

o o 

0 

F° — 

c? 

o 

1 

o 

CD 

"CiO 

1 

, r — 

1 

O 



( 2 . 20 ) 


The unassembled closed-loop system is an optimal control system since each 
subcontroller is optimal for its associated substructure. 

The last step is to assemble the subcontrollers by using the same coupling 
scheme as that used for assembling the substructures. The assembled controller 



10 


for the assembled structure is represented by 

Sq = (A + BG ® - F*C)q + F*y 
u = G®q 

with 

F® = T t F°, G ® = G°T 


( 2 . 21 ) 


( 2 . 22 ) 


where superscript ® denotes that the controller is not optimal but is considered 
as suboptimal. 

One advantage of using Substructural Controller Synthesis to design a 
controller is that an SCS controller is highly adaptable. For a structure with 
varying configuration or varying mass and stiffness properties, like some space 
structures, the Substructural Controller Synthesis method may be especially 
efficient. The SCS controller can be updated economically by simply carrying 
out redesign of subcontrollers associated with those substructures that have 
changed. On the other hand, for a controller based on a centralized design 
scheme, a slight change of the structure may require a full-scale controller 
redesign. 

An LQGSCS Algorithm is presented in the paper to summarize the Sub- 
structural Controller Synthesis procedure. Also a two-component planar truss 
structure with non-colocated actuators and sensors is used to demonstrate the 
applicability of the proposed method. 


2.3 A Review of Model Reduction Methods for Struc- 
tural Control Design (Ref. 10) 

In this paper, several frequently used model reduction methods are briefly 
reviewed. The methods reviewed include: modal truncation, balanced model 
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redaction, balanced gain approach, Krylov model reduction, and Ritz vectors 
and mixed- mode method. 

Among a myriad of existing model reduction methods for structural dy- 
namics systems, modal truncation may be the most popular approach. A modal 
representation has many advantages: modal frequencies represent resonances 
of the structure, equations of motion are uncoupled implying a saving of com- 
putation time, and modal data can be identified and validated by vibration 
test. However, selection of modes to be retained in the reduced model may not 
be an easy task. The simplest approach would be to include all modes within 
the frequency range of interest. For a large space structure with closely- spaced 
frequencies, this simplest approach may produce a reduced model whose size 
is still too large to handle. 

An efficient modal truncation criterion is based on balanced singular val- 
ues. If the structure frequencies are sufficiently separated and modal damping 
is very small, then modal representation of a structural dynamics system is 
approximately balanced. The approximate balanced singular values are calcu- 
lated by 


where p, is the z'-th row of force distribution matrix, u t - and ty, are the z’-th 
columns of the displacement and velocity sensor distribution matrices, and 
A, and 77, are the z-th modal frequency and damping ratio. Balanced model 
reduction is performed on modal coordinates. 

Other than balanced singular values, balanced gains also can serve as a 


( 2 . 23 ) 
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basis for modal truncation. Balanced gains are defined as 

9i = 

The balanced gain approach is an optimal modal truncation method in the L 2 
sense. It is closely related to the modal cost analysis. 

In addition to normal modes, there are other Ritz vector superposition 
methods for dynamic analysis of structures. In Ref. [4], Su and Craig presented 
a Ritz vector method called Krylov model reduction. Krylov vectors are system 
static modes generated by a recurrence procedure. For undamped structural 
dynamics systems, the Krylov procedure is 

Q j+ 1 = K-'MQj (2.25) 



For damped systems, the Krylov procedure is 

/ QU \ - f ~ K ~' D -K~'M 

1 “ l I 0 



(2.26) 


Krylov reduced models match a set of system parameters called low-frequency 
moments. 


Other than using only normal modes or only Ritz vectors, a mixed-mode 
method combines some dominant normal modes and some static modes in the 
basis for model reduction. Static modes are a system’s deflection shapes associ- 
ated with imposed force distribution vectors. Krylov vectors can be considered 
to be system static modes. Recently, several numerical experiments have shown 
that by augmenting a modal basis with some Ritz vectors, fidelity of the re- 
duced model can be substantially improved [13]. 
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In this paper, a plane truss structure with closely-spaced frequencies and 
light modal damping is used to compare different reduced-order models. Open- 
loop comparison includes the L 2 error norm of the impulse response function 
and approximation of the output frequency response function. The control de- 
sign comparison includes closed-loop stability and control performance. For the 
example studied, the balanced gain approach produces reduced-order models 
that approximate the impulse response better than any other truncation cri- 
terion. However, for closed-loop comparison, modal truncation by preserving 
the lowest frequency modes yields more stable closed-loop designs than other 
methods. 

2.4 Recent Literature on Structural Modeling, Identi- 
fication, and Analysis (Ref. 12) 

This paper, which surveys literature from 1980 to the present related to 
the topics of modeling, identification, and analysis of large space structures, 
provides a list of over 240 references. The topic of “Mathematical Modeling of 
Large Space Structures” includes sections dealing with the following subjects: 

• Continuum Models 

• Model Order Reduction 

• Substructuring; Component Synthesis 

• Computational Techniques 


• Waves vs Modes 
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• Localization: Local vs Global Behavior 

• Nonlinearity: Joints 

• Damping 

• Modeling Errors and Parameter Uncertainty 

Under the topic “System Identification; Model Verification and Model 
Updating” the following subjects are surveyed: 

• Experimental Modal Analysis 

• Model Verification: Model Error Localization and Model Updating 

• Scale Modeling of Large Space Structures 

• Damage Detection 


• On-Orbit System Identification 



Chapter 3 


PROJECTED FUTURE WORK 


During the second year of this project, the research will continue to focus 
on investigation of important issues in structure/control interaction problems, 
and will specifically focus on the following topics: 

1. To address stability issues in substructure-based control design, and to 
develop an augmented global controller in order to reduce the influence of 
the substructures on each other in the closed-loop system, such that sta- 
bility and performance of the substructure-based controllers is enhanced. 

2. To develop a reduced-order observer in second-order form, which is to be 
used as the feedback basis in the control design for large scale structural 
dynamics systems. 

3. To incorporate the Krylov model reduction and component mode syn- 
thesis method with the substructural control design approach, and to 
illustrate the feasibility and efficiency of the Substructural Controller 
Synthesis method when it is used to design controllers for practical large 
flexible space structures. 

4. To apply the Krylov model reduction and the substructural controller 
synthesis method to the controller design for multi-flexible-bodies and 
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articulated mechanical systems, which involve reduced-order modeling, 
precision-pointing, and vibration suppression. 
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